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The recently reported measurements of the CP asymmetry a^K by the BABAR 
and BELLE collaborations, obtained from the rate differences in the decays B° — > 
(J/tpKs), (J/ipKi,) etc., and their charge conjugates, are in good agreement with the 
standard model (SM) prediction of the same, resulting from the unitarity of the CKM 
matrix. The so-called minimal flavour violating (MFV) supersymmetric extensions 
of the standard model, in which the CKM matrix remains the only flavour changing 
structure, predict cl^k similar to the one in the SM. With the anticipated precision in 



; a^K an d other CP asymmetries at the B factories and hadron colliders, one hopes to 

■ pin down any possible deviation from the SM. We discuss an extension of the MFV- 

super symmetric models which comfortably accommodates the current measurements 
of the CP asymmetry a^jc, but differs from the SM and the MFV-supersymmetric 
models due to an additional flavour changing structure beyond the CKM matrix. We 



suggest specific tests in forthcoming experiments in B physics. In addition to the CP- 
asymmetries in .B-meson decays, such as cl^k and a n7T , and the mass difference AM S in 
the B® - B® system, we emphasize measurements of the radiative transition b — ► d'y as 
sensitive probes of the postulated flavour changing structure. This is quantified in terms 
of the ratio R(py/K*j) = 2B(B° -> p°-/)/B{B° -> K*°i), the isospin violating ratio 
A ±0 = B(B^ — ► p ziz "f)/2B(B° — > p 0, y) — 1, and the CP-asymmetry in the decay rates 
for B + — > p + 7 and its charge conjugate. Interestingly, the CKM-unitarity analysis in 
the Extended-MFV model also allows solutions p < for the Wolfenstein parameter, 
as opposed to the SM and the MFV-supersymmetric models for which only p > 
solutions are now admissible, implying 7 > tt/2, where 7 = — argV^. Such large 
values of 7 are hinted by the current measurements of the branching ratios for the 
decays B — > tttt and B — > Ktt. 
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1 Introduction 



With the advent of the B-factory era, the principal focus in flavour physics is now on mea- 
suring CP-violating asymmetries, which will determine the inner angles a, (3, and 7 of the 
unitarity triangle (UT) in the Cabibbo-Kobayashi-Maskawa (CKM) theory Q. A begin- 
ning along this road has already been made through the impressive measurements of sin 2/? 
by the B-factory experiments BABAR [0 and BELLE ||, following earlier leads from the 
OPAL H], CDF |5|,|[J, and ALEPH [0] collaborations. The principal decay modes used in 
the measurement of sin 2/3 are B® — > J/ipKs, B% — > ip{2S)Ks, B% — > J/i/jKl, and their 
charge conjugates. Concentrating on the decays B%/B® — > J/ipK s , the time-dependent 
CP-asymmetry a^,K- s (t) can be expressed as follows: 



a>ipK s if) 



r(B°(t) -> j/^) - r(gg(t) -» J/y^) 
r(B°(t) - + r(BS(t) - j/^k s ) 

A$ P cos(AM Bd t) + ^ sin( AM Bd t) , (1) 



where the states B%(t) and are understood as evolving from the corresponding initial 

flavour eigenstates (i.e., at t — 0), and AM^ is the mass difference between the two mass 
eigenstates of the B® -B% system, known very precisely, thanks in part due to the BABAR || 
and BELLE || measurements, and the present world average is AM^ d = 0.484 ± 0.010 
(ps)" 1 The quantities A$p and .4^ x are called the direct (i.e., emanating from the 

decays) and mixing-induced CP-asymmetries, respectively. Of these, the former is CKM- 
suppressed - a result which holds in the SM. The expectation Aq P /Aqp <C 1 is supported 
by present data on direct CP- asymmetry in charged B-decays, B^ — > J/i/jK^, yielding an 
upper bound on Aq P which is already quite stringent @,|H- Hence, we shall assume that 
direct CP-asymmetry in a^- s (t) is negligible and neglect the first term on the r.h.s. of 
Eq. Recalling that ^4™p x * s a pure phase, one has in the SM Aq P x = r lJ/i>K s O"ipK s , with 
Tju/^Ks) = ~ 1 bdng the intrinsic CP-parity of the J/ipKs state, Eq. flU) simplifies to 

a^K s {t) = -cifK s sin(AM Bd t) ■ (2) 

This relation is essentially free of hadronic uncertainties. Hence, a measurement of the left- 
hand-side allows to extract a^K s cleanly. Note that in the SM a^K s = sin 2(3. The present 
world average of this quantity is 0-[7|] 

a^Kg = 0-79 ±0.12 , (3) 

which is dominated by the BABAR {a^ Kg = 0.59 ± 0.14(stat) ± 0.05(syst)) and BELLE 
{a>i>K s = 0.99 ± 0.14(stat) ± 0.06(syst)) || results. We note that the current world average 
(which includes a scale factor following the Particle Data Group prescription [[□]]) based 
on all five experiments yields a value of a^K s which is different from a null result by more 
than six standard deviations. To test the consistency of the SM, the current experimental 
value of a^K s i n Eq. i^j is to be compared with the indirect theoretical estimates of the same 
obtained from the unitarity of the CKM matrix. These latter values lie typically in the range 



a 



SM 

ipK s 



0.6 — 0.8 (at 68% C.L.) |12|-p!9[], where the spread reflects both the uncertainties in 



1 



the input parameters and treatment of errors, with most analyses yielding (iTk s — 0.70 as 
the central value of the CKM fits. We conclude that he current measurements of 

CL-tpKs <^re 

in good agreement with its indirect estimates in the SM. 

The consistency of the SM with experiments on CP-violation in B-decays will come 
under minute scrutiny, with greatly improved accuracy on a^K s and measurements of the 
other two angles of the UT, a and 7 at the e + e~ and hadronic U-factories. In addition, a 
large number of direct CP-asymmetries in charged and neutral -B-decays, as well as flavour- 
changing-neutral-current (FCNC) transitions in B- and f^-decays, which will be measured 
in the course of the next several years, will greatly help in pinning down the underlying 
theory of flavour physics. It is conceivable that precision experiments in flavour physics may 
force us to revise the SM framework by admitting new interactions, including the possibility 
of having new CP-violating phases. Some alternatives yielding a lower value of a^K s than 
in the SM have already been entertained in the literature |20|-p2|. With the experimental 



situation now crystallized in Eq. (^), it now appears that the CP-asymmetry cl^k s has a 
dominant ly SM origin. 

In popular extensions of the SM, such as the minimal supersymmetric standard model 
(MSSM), one anticipates supersymmetric contributions to FCNC processes, in particular 
AMB d , AMb s (the mass difference in the -B® system), and €k, characterizing Aq™ 
in the K° -K° system. However, if the CKM matrix remains effectively the only flavour 
changing (FC) structure, which is the case if the quark and squark mass matrices can be 
simultaneously diagonalized (equivalently, the off-diagonal squark mass matrix elements are 
small at low energy scale), and all other FC interactions are associated with rather high 
scales, then all hadronic flavour transitions can be interpreted in terms of the same unitarity 
triangles which one encounters in the SM. In particular, in these theories a^K s measures 
the same quantity sin 2(3 as in the SM. These models are usually called the minimal flavour 
violating (MFV) models, following Ref. p3| . Despite the intrinsic dependence of the mass 



differences AM^, AMb s , and €k on the underlying supersymmetric parameters, the MFV 
models remain very predictive and hence they have received a lot of theoretical attention 
lately [l5|,^3|-|2j|. To summarize, in these models the SUSY contributions to AMe d , AMb s , 



and 6k have the same CKM-dependence as the SM top quark contributions in the box 
diagrams (denoted below by C^ tt ). Moreover, supersymmetric effects are highly correlated 
and their contributions in the quantities relevant for the UT-analysis can be effectively 
incorporated in terms of a single common parameter / by the following replacement 0, : 

e K , AM Bs AM Bd , a^ Ks : C?* - Cf"(l + /) . (4) 

The parameter / is positive definite and real, implying that there are no new phases in 
any of the quantities specified above. The size of / depends on the parameters of the 



supersymmetric models and the model itself |3(J-[33|]. Given a value of /, the CKM unitarity 



fits can be performed in these models much the same way as they are done for the SM. 
Qualitatively, the CKM-fits in MFV models yield the following pattern for the three inner 
angles of the UT: 

^MFV^SM. 7 MFV <7 SM. ^MFV > Q SM (5) 

For example, a recent CKM- fit along these lines yields the following central values for the 
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three angles 
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/ = (SM) : («,/?, 7 ) central = (95°, 22°, 63°) , 

/ = 0.4 (MFV): (a,/3,7)cent r ai = (112 o ,20 o ,48 o ) . (6) 

leading to (sin 2/?)^ tral ~ 0.70 and (sin 2/3)^^ ^ .64. Thus, what concerns sin 2/5, the 
SM and the MFV models give similar results from the UT-fits, unless much larger values for 



the parameter / are admitted which, as argued in Refs. ||30|-p3l and in this paper, is unlikely 
due to the existing constraints on the MFV-SUSY parameters. 

However, in a general extension of the SM, one expects that all the quantities appearing 
on the l.h.s. in Eq. ([|) will receive independent additional contributions. In this case, the 
magnitude and the phase of the off-diagonal elements in the B® -B® and B° s - B® mass 
matrices can be parametrized as follows [34, 35]: 



M l2 (B d ) = {BdlH 2 ^ lBd) = &»<Mg<(B d ) , 

M 12 (B S ) = ™ = ry^M™(B s ) . (7) 



where r d (r s ) and 9 d ( 6 d ) characterize, respectively, the magnitude and the phase of the new 
physics contribution to the mass difference AM^ (AMgJ. It follows that a measurement 
of a^K a would not measure sin 2(3, but rather a combination sin2(/5 + 9 d ). Likewise, a 
measurement of the CP asymmetry in the decays B^ — > tttt and its charge conjugate, a 7Tn , 
(assuming that the penguin contributions are known) would not measure sin 2a, but rather 
sin2(a— 9 d ). Very much along the same lines, the decay B s — > J / ip<p and its charge conjugate 
would yield a CP asymmetry ~ — sin(<5 — 9 S ), where 5 ~ 1/2A 2 ~ 0.02 in the SM, and 



A ~ 0.22 is one of the four CKM parameters in the Wolfenstein representation | 36fl . Thus, 
the phase 9 S could enhance the CP-asymmetry bringing it within reach of the LHC- 
experiments [j3~7fl . In this scenario, one also expects new contributions in Mi 2 {K°), bringing 



in their wake additional parameters (r e , 9 e ). They will alter the profile of CP- violation in the 
decays of the neutral kaons. In fact, sizable contributions from the supersymmetric sector 
have been entertained in the literature, though it appears now unlikely that ex and/or e' K /ex 
(which is a measure of direct CP violation in the neutral Kaon decays) are saturated by 



supersymmetry |p8| , |39|| . 



It is obvious that in such a general theoretical scenario, which introduces six a priori 
independent parameters, the predictive power vested in the CKM-UT analysis is lost. We 
would like to retain this predictivity, at least partially, and entertain a theoretical scenario 
which accommodates the current measurement in flavour physics, including the recent mea- 
surements of a^K S i but admits additional flavour structure. A model which incorporates 
these features is introduced and discussed in section 2, using the language of minimal inser- 
tion approximation (MIA) ||0| in a supersymmetric context. In this framework, gluinos are 
assumed heavy and hence have no measurable consequences for low energy phenomenology. 
All FC transitions which are not generated by the CKM mixing matrix are proportional to 
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the properly normalized off-diagonal elements of the squark mass matrices: 



, - , .l> ( M ij)AB , , 



where i,j = 1, 2, 3 and A,B = L, R. We give arguments why we expect that the dominant 
effect of the non-CKM structure contained in the MIA-parameters is expected to influence 
mainly the b —>■ d and s — > d transitions while the b — > s transition is governed by the 
MFV-SUSY and the SM contributions alone. For what concerns the quantities entering in 
the UT analysis, the following pattern for the supersymmetric contributions emerges in this 
model: 

AM Bs : Cf « -> Cf "(1 + f) (9) 
e K , AM Bd , a+ Ka : Cf* - C?«{1 + /) + Cf r = C?« (1+f + g) (10) 

where the parameters / and g = gn + igi represent normalized (w.r.t the SM top quark Wtt) 
contributions from the MFV and MIA sectors, respectively. Thus, in the UT-analysis the 
contribution from the supersymmetric sector can be parametrized by two real parameters / 
and gn and a parameter gi, generating a phase 9^, which is in general non-zero due to the 
complex nature of the appropriate mass insertion parameter. We constrain these parame- 
ters, taking into account all direct and indirect bounds on the supersymmetric parameters, 
including the measured rates for b — > 57 decay |41|-|43|, (g — 2) M from the Brookhaven exper- 
iment [03] , and the present bound on the b — > d^y transition, following from the experimental 



bound on the ratio of the branching ratios R^fry/K*^) = 2B(B° — > p°'y)/B(B — > K*°j) Q . 
We do not include the quantity e' K /eK in our analysis, despite its impeccable measure- 



ment by the NA48 J|6| and KTeV |47| Collaborations, yielding the present world average 
Kee' K /eK = (1.7 ± 0.2) x 10 -3 , due to the inherent non-perturbative uncertainties which 
have greatly reduced the impact of the e' K /eK measurement on the CKM phenomenology 
(see, for a recent review, Ref. fTi 



This model, called henceforth the Extended-MFV model, leads to a number of testable 
consequences, some of which are common with the more general scenarios discussed earlier 
in the context of Eq. (|7|) []5^|55|. Thus, for a certain range of the argument of the MIA 



parameter, this model yields a^K s < a ^,^ MFV - F° r other choices of the model parameters, 
this model yields a higher value for this CP asymmetry. A precise measurement of a^,K s 
would fix this argument (= 9a) and we show its allowed range suggested by the current data. 
Likewise, the CP-asymmetry a n7T = sin 2(a— 9 4) will be shifted from its SM- value, determined 
by 9d- All b — > d transitions (leading to the decays such as b — >• dj, b — > d£ + £~, — > 
where £ ± = e^, /i 1 * 1 , t^, and the ratio of the mass differences AM^/AM^J may turn out to 
be significantly different from their SM and MFV counterparts. To illustrate this, we work 
out in detail the implications for the exclusive decays B — > p°7 and B ± — ► p ± 7, concentrating 
on the (theoretically more reliable) ratios R(pj/K*j) = 2B(B° — > p°'y)/B(B° — > K*j), the 
isospin violating ratio A ±0 = B{B ± p ±/ y) /2B(B° — ► p°j) — 1, and direct CP-asymmetry in 
the decay rates for B — > p~7 and its charge conjugate. We also find that the fits of the CKM 
unitarity triangle in the extended MFV model, characterized by Eqs. (|9[) and ( ]TXJ| ) above, 
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admit both p > and p < solutions, where p is one of the Wolfenstein parameters |3~6[|. We 



illustrate this by working out the predicted values of 7 (and a) and AM# s in this model for 
some specific choice of the parameters. This is in contrast with the corresponding fits in the 
SM and the MFV-MSSM models, which currently yield p > at 2 standard deviations in 
the SM, with the significance increasing in the MFV-MSSM models. The allowed CKM-fits 
in the extended-MFV model with p < imply in turn 7 > tt/2. We note that such large 



values of 7 are hinted by phenomenological analyses [58, 49[ of the current measurements 
of the branching ratios for the decays B — > Kit, tttt [^0|-|52]|. However, this inference is 



not yet convincing due to the present precision of data and lack of a reliable estimate of 
non-perturbative final state interactions in these decays. 

This paper is organized as follows: In section 2, we give the outline of our extended-MFV 
model. The supersymmetric contributions to the quantities of interest (ex, AMb s , AM^, 
a*l>K s ) and -R(p7 / 'K*y) are discussed in section 3, where we also discuss the impact of the 
(g — 2)^ experiment on our analysis. Numerical analysis of the parameters (/, \g\), taking into 
account the experimental constraints from the b — > 37, b — > dj and (g — 2)^, is presented in 
section 4. A comparative analysis of the unitarity triangle in the SM, MFV and the Extended- 
MFV models is described in section 5, where we also show the resulting constraints on the 
parameters and gi and the CP-asymmetry a$K s - The impact of the Extended-MFV 
model on the b — > d'-f transitions are worked out in section 6. Section 7 contains a summary 
and some concluding remarks. The explicit stop and chargino mass matrices are displayed 
in Appendix A and some loop functions encountered in the supersymmetric contributions 
are given in Appendix B. 



2 Outline of the model 

The supersymmetric model that we consider is a generalization of the model proposed in 
Ref. |[23|| , based on the assumptions of Minimal Flavour Violation with heavy squarks (of 
the first two generations) and gluinos. The charged Higgs and the lightest chargino and 
stop masses are required to be heavier than 100 GeV in order to satisfy the lower bounds 
from direct searches. The rest of the SUSY spectrum is assumed to be almost degenerate 
and heavier than 1 TeV. In this framework the lightest stop is almost right-handed and the 
stop mixing angle (which parameterizes the amount of the left-handed stop present in the 
lighter mass eigenstate) turns out to be of order 0(Mw/Mg) ~ 10%; for definiteness we will 
take |0 f | < vr/10. 

The assumption of a heavy (> 1 TeV) gluino totally suppresses any possible gluino- 
mediated SUSY contribution to low energy observables. On the other hand, the presence 
of only a single light squark mass eigenstate (out of twelve) has strong consequences due to 
the rich flavour structure which emerges from the squark mass matrices. As discussed in the 
preceding section, adopting the MIA- framework ||40| , all the FC effects which are not gen- 
erated by the CKM mixing matrix are proportional to the properly normalized off-diagonal 
elements of the squark mass matrices (see Eq. (|§D). In order to take into account the effect 
of a light stop, we exactly diagonalize the 2x2 stop system and adopt the slightly different 
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MIA implementation proposed in Ref. ||53|| . In this approach, a diagram can contribute siz- 
ably only if the inserted mass insertions involve the light stop. All other diagrams require 
necessarily a loop with at least two heavy (> 1 TeV) squarks and are therefore automatically 
suppressed. This leaves us with only two unsuppressed flavour changing sources other than 
the CKM matrix, namely the mixings ul — ^2 (denoted by $u L t 2 ) an d c_l — ^2 (denoted by 
bc L i 2 )- We note that 5 fi r 2 and 6 5 f 2 are mass insertions extracted from the up-squarks mass 
matrix after the diagonalization of the stop system and are therefore linear combinations of 
( s is)lr: ( S is)ll and of OfeOLj, ( S 23,)ll, respectively. 

Finally, a comment on the normalization that we adopt for the mass insertions is in 
order. In Ref. [Qj it has been pointed out that ((^3)^ must satisfy an upper bound of order 



2m t /Mq- in order to avoid charge and colour breaking minima and directions unbounded 
from below in the scalar potential. We normalize the insertions relevant to our discussion so 
that, in the limit of light stop, they automatically satisfy this constraint: 



M-2 ? \V+ 

r _ u(c) L t 2 I y tays) \ /-.^.n 

1V1 t 2 1V1 q V td(s) 



This definition includes the phase of the CKM element V t d( s )- m this wa Y; deviations from 
the SM predictions, for what concerns CP violating observables, will be mainly associated 
with complex values of the mass insertion parameters. For instance, as we will argue in 
the following, the CP asymmetry a^x s in the decay B — > J/ipK s can differ from the SM 
expectation only if axgS^^ 7^ 0. In general, the two phases must not be aligned with the 
respective SM-phases entering in the box diagram: 

-*7ibt (i2 » 

The insertion $c L t 2 characterizes the b —* s transitions and it enters in the determination of 
the B s — B s mass difference, the b — > 57 decay rate, and observables related to other FCNC 



decays such as b — > s£ + £~~. For what concerns the b — > S7 decay, previous analyses [po 
pointed out that contributions proportional to this insertion can be as large as the SM one. 
The experimental results for the inclusive branching fraction are 

(3.19 ± 0A3 atat ± 0.27 syst ) x 10" 4 CLEO [41] 

B(B -> X sl ) = { (3.11 ± 0.80 stat ± 0.72 3 , st ) x 10" 4 ALEPH || (13) 



(3.36 ± 0.53 stat ± A2 ayst ± (° f 4 ) model ) x 10" 4 BELLE @ 



Combining these results and adding the errors in quadrature we obtain the following world 
average for the inclusive branching ratio 

B(B -> X sl ) = (3.22 ± 0.40) x 10~ 4 (14) 

yielding the following 95% C.L. experimentally allowed range 

2.41 x 10" 4 < B(B -> X s7 ) < 4.02 x 10" 4 . (15) 
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Using the LO theoretical expression for this branching ratio, the following bound is obtained: 



0.30 < 



C 7 ff ' LO {m b ) <0.40 (16) 



where L (m b ) is the relevant Wilson coefficient evaluated in the LL approximation and 
has the value -0.316 in the SM. Analyses of the NLO SM [§§ and SUSY [T2|,|S7}fU[ con- 
tributions (for the latter only a limited class of SUSY models were considered) showed that 
the LO result can receive substantial corrections. Notice that the NLO analysis presented 
in Refs. p3 , [57]HH! is valid exactly for the class of models that we consider over all the 



SUSY parameter space (including the large tan/3 region). Implementing their formulae and 
allowing the SUSY input parameters to vary within the range that we will further discuss 
in section f|, we find that, up to few percents, the branching ratio is given by 

B{B^X sl ) = 22.23 \C^ ff ' NLO (m b ) - 0.06l] 2 + 0.264, (17) 
Cj ff ' NLO (m b ) = -0.175 + 0.666 C^ f ' NLO (M w ) + 0.093 Cl ff ' NLO (M w ) . (18) 
The explicit expressions for C^/ ,NLO (M w ) can be found in Refs. P3],|57|-pH. Combining 



Eq. (\V7\ j and the bound (|15|), we obtain 

- 0.35 < C e 7 ff ' NLO (m h ) < -0.24 or 0.36 < C e 7 ff ' NLO (m b ) < 0.49 . (19) 



Notice that in Refs. |58],|5£| it is pointed out that in order to get a result stable against 
variations of the heavy SUSY particles scale, it is necessary to properly take into account all 
possible SUSY contributions and to resum all the large logarithms that arise. The inclusion 
of the insertion d> 5 r in this picture, in particular, should not be limited to the LO matching 
conditions but should instead extend to the complete NLO analysis. This program is clearly 
beyond the scope of the present paper. Moreover, one finds that, including the NLO correc- 
tions, the SM almost saturates the experimental branching ratio. In view of this we choose 
not to consider S dL ^ 2 in our analysis. Of course, the SUSY contribution from the MFV sector 
is still there, but it is real relative to the SM. The assumption of neglecting 5 5l ^ 2 will be 
tested in CP-asymmetries Acp(b ~~ > s l) an d Acp(B — > K*j) at the B-factories. Notice that 
the exclusion of 8 5 r from our analysis introduces strong correlations between the physics 
that governs b — > d and b — > s transitions, such as the ratio AM B J AM Bd , which would 
deviate from its SM (and MFV model) values. 

The free parameters of the model are the common mass of the heavy squarks and gluino 
(Mg), the mass of the lightest stop (Mf 2 ), the stop mixing angle (%), the ratio of the two 
Higgs vevs (tan/3,5 Q), the two parameters of the chargino mass matrix (/i and M 2 ), the 
charged Higgs mass (M H ±) and &u L ii- All these parameters are assumed to be real with the 
only exception of the mass insertion whose phase in not restricted a priori. In this way we 
avoid any possible problem with too large contributions to flavour conserving CP violating 
observables like the electric dipole moments of the leptons, hadrons and atoms. 

In the next section we analyze the structure of the SUSY contributions to the observables 
related to the determination of the unitarity triangle, namely e#, AM djS and a^K s - 

^We adopt the notation /3g in order not to generate confusion with the inner angle of the unitarity triangle 
which is denoted by (3 
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3 SUSY contributions 



The effective Hamiltonian that describes AF = 2 transitions can be written as 



n 



AF=2 
eff 



C 2 /If 2 



(2tt) 2 y ' iqi tq2 
+ CM q% L qf R ■ 4lQ.xr) + h - c - 



(20) 

where a, (3 are colour indices and (qi, q 2 ) = (s, d), (6, d), (b, s) for the K, B d and B s systems 
respectively. 

In this framework, as previously explained, gluino contributions are negligible; therefore, 
we have to deal only with charged Higgs (which obeys the SM CKM structure) and chargino 
mediated box diagrams. Let us comment on the latter. The dominant graphs must involve 
exclusively the lightest stop eigenstate since the presence of an heavy (> 1 TeV) sparticle 
would definitely suppress their contribution. Moreover, Feynman diagrams that contribute 
to C 3 are substantially suppressed with respect to diagrams that contributes to C\. In fact, 
for what concerns the B system, the vertices 5r — H — t 2 and — H — i 2 are proportional, 
respectively, to m b sin 9^ / (V2 M w cos f3s) and m t cos9i/(V2M w smf3 s )- Their ratio is thus 
of order (mb/rrit tan fis tan 9i) 2 which, even in the large tan /3s regime, is damped and not 
much larger than O(0. 1). In the K system m& must be replaced by m s and the suppression 
is even stronger. Notice that in frameworks in which the split between the two stop mass 
eigenstates is not so marked, this argument fails. Diagrams mediated by the exchange of 
both stops must be considered and it is possible to find regions of the parameter space (for 
large tan (3s) in which SUSY contributions to C3 are indeed dominant [3S|. We note that, 
the tan 4 (3$ enhanced neutral Higgs contributions to the coefficients 6*2,3, whose presence is 
pointed out in Ref. U60| , do not impact significantly for the range of SUSY parameters that 
we consider (tan (3s < 35 and \9^\ < it/ 10). 

The total contribution to C\ can thus be written as 



C[ ot {M w ) = CY{M W ) + C^(M W ) + C*{M W ) + Cf n {M w ) 



MIi 



The explicit expressions for the various terms are: 



Cf = 



1 ~ 4 tan 4 /3 5 

-7 — TIT Y 2\l,%H,Xt,Xt 

tan z /3s 



(V V* \ 2 V V* 

1 \v tqi v t *J 1 v tqi v t * 2 1 

i~y> 2 / >-» 2 



(21) 



(22) 



Yi{xh, x h , x t , x t ) + - \ - Yi(l, x H , x t , x t ) 
2 taxr (3$ 



c? = E 



~ m t cos 9~ t . ~ 
V2M W sm (j s 



i cos 9i 



x ^ / i( a; t? x t, x Xi , x Xj/ 



m t 

j2 y/2M w ahiPs 



-sin 0i% 



(23) 



(24) 
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MI 



ud 



Id 




m t cos 61 



y/2M w sin (3 S 

X ^l { x q> x h X Xti X Xj) ^u L h 
y^MI X 2 
°1 °u L h 



mnOfe) [V j2 



m t cos 61 
V2M W sin (3 S 



(25) 



where C 1 a ^ = G(x a , xp) and x a = m^/M^. The contributions to and Cf 11 come from 
the Feynman diagrams shown in Fig. |I|. The conventions we adopt for the chargino mass 
matrix, the exact definition of the stop mixing angle and the explicit expressions for the 
loop functions G, Y%, Y% and Y± can be found in the appendices. Notice that, if we restrict 
to the Bd and B s cases, one may neglect the terms C± tc and C^ cc in Eq. (^) as they are 
suppressed by small CKM factors. In the K system, on the other hand, it is necessary to 
consider all the terms. 

Eq. (|25|) describes the impact of a non-zero mass insertion on the Bd and K systems. 
The corresponding contribution to the B s system is obtained via the substitution 



(26) 



This implies that the impact of this diagram on the B s system is reduced by a factor 
|KisVjd/Vj s | 2 ~ 0.0016 with respect to the Bd and K systems. Since, as we will show in 
the forthcoming analysis, Cf n is not likely to exceed by more than twice the SM contri- 
bution, it is clear that any effect in the B s system from the mass insertion is completely 
negligible. 

As already discussed in the introduction, the following structure of the SUSY contribu- 
tions emerges in the class of model described above: 



V ud 


2 


v us 


Vtd 




V ts 



AM Bs 
AMe d , a^Ks 



(jWtt 
(jWtt 



cr*(i + /) 

Cf tt (l + /) + Cf J = Cf * (1 + / + 9) 



where the parameters / and g represent normalized contributions from the MFV and MIA 
sectors, respectively, 



/ 
9 



+ cf)/cl 

gR + igi = C 1 



Wtt 



S 2 r 
ULt2 



Wtt 



(27) 
(28) 



The impact of the SUSY models on the observables we are interested in is then parametrized 
by three real parameters /, gn and gj. We recall here that in the limit g — > 0, the above 

23.1213] for the minimal flavour violation 



parametrization reduces to the one given in Ref. 
and the MSSM cases (More generally, for all models in which the CKM matrix is the only 
flavour-changing structure). The absence of any CKM phase in Eq. (^) as well as in the 
definition of g reflects the definition of the mass insertion given in Eq. (pT|). 

Note that / and g are functions of the SUSY parameters that enter in the computation of 
many other observables that are not directly related to CP violation. This implies that it is 
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possible to look for processes that depend on the same SUSY inputs / and g. For example, 
the presence of non trivial experimental bounds on the \AB = 2|, AQ = transitions can 
induce interesting correlations with the UT analysis. Likewise, the inclusive radiative decay 
B — ► Xsj and (g — 2) M , the anomalous magnetic moment of the muon, are susceptible to 
supersymmetric contributions. On the other hand, it is necessary to search for observables 
that can put constraints on the insertion 5 &L f 2 . In this context, the transitions b —>■ dj, 
B° — > p°7 and B ± — > p ±r y are obvious places to look for 8u L t 2 -related effects. 

Concerning now the (g — 2) M constraint, we recall that the Brookhaven Muon (g — 
2)-collaboration has recently measured with improved precision the anomalous magnetic 



moment of the positive muon. The present world average for this quantity is [13] 

a M+ (exp) = 116592023(15) x 10~ 10 . (29) 

The contribution to a M + in the SM arises from the QED and electroweak corrections and 
from the hadronic contribution which includes both the vacuum polarization and light by 
light scattering |H]] . The error in the SM estimate is dominated by the hadronic contribution 
to and is obtained from a(e + e~ — > hadrons) via a dispersion relation and perturbative 
QCD. The light-by-light hadronic contribution is, however, completely theory-driven. Sev- 
eral competing estimates of a^+ d exist in the literature, reviewed recently in Ref. [|62[| . We 
briefly discuss a couple of representative estimates here. 

In order to minimize the experimental errors in the low-s region, Davier and Hocker 
supplemented the e + e~ — ► 7r + 7r~ cross section using data from tau decays. Using isospin 
symmetry they estimate p3| 



a h t = (692.4 ± 6.2) x 10~ 10 . (30) 



An updated value of a^+ d using earlier estimates of Eidelman and Jegerlehner supple- 
mented by the recent data from CMD and BES collaborations yields [^5| 

al% d = (698.75 ± 11.11) x 10~ 10 . (31) 

Alternatively, calculating the Adler function from e + e~ data and perturbative QCD for the 
tail (above 11 GeV) and calculating the shift in the electromagnetic fine structure constant 
Aa had in the Euclidean region, Jegerlehner quotes p5| 



a Jf = (697.4 ± 10.45) x 10~ 10 . (32) 

These estimates are quite compatible with each other, though the errors in Eqs. (|3~ID and (|32j) 
are larger than in Eq. (|30"D . Adopting the theoretical estimate from Davier and Hocker [|63|] . 
one gets 

a M+ (SM) = 11659163(7) x 10~ 10 , (33) 

yielding 

5a^+ = a M+ (exp) - a M +(SM) = +43(16) x 10" 10 , (34) 

which is a 2.6 a deviation from the SM. Using, however, the estimate from Jegerlehner in 
Eq. (H), gives 

Sa^ = +37(18) x 10~ 10 , (35) 



10 



which amounts to about 2 a deviation from the SM. Thus, on the face-value, there exists a 

2) M and SM-estimates. 

ii nnrl a 



2 to 2.6 a discrepancy between the current experiments on (g 



In SUSY theories, 



receives contributions via vertex diagrams with x ft an d X 



loops pl|,p6|-|7t: 



The chargino diagram strongly dominates over almost all the parameter 
space. The chargino contribution is |66| (see also Ref. |67] for a discussion on CP violating 
phases): 



92 



m: 



2 2 



8tt 2 Ml 



E 



Ml Re {Uf 



n) 



2M2 /C os 2 /3 s 



+ ^ (Vl) 



M Xi Re(U l2 V a ) 
V2M W cos (3 



[X; 



(36) 



is 



where Xi = M^./M~ and the loop functions Fi and F 3 are given in the appendices. Eq. ( p6| 

dominated by the last term in curly brackets whose sign is determined by sign[Re([/i2Vii)] = 
— sign[Re(/i)] (note that we have M Xl < M X2 ). Taking into account that the Brookhaven 



experiment implies <5a M + > at 2 to 2.6 a, it is clear that the fi > region is currently 
favoured. 



Finally let us focus on & - 
limit on the inclusive decay b - 
exclusive branching ratios are 



g?7 decays. To the best of our knowledge, there is no direct 
-> d'y. The present experimental upper limits on some of the 



B{B° -> p° 7 ) < 0.56 x 1(T 5 (90% C.L.) [45] , 
B(B + -> p+ 7 ) < 1.3 x 10- 5 (90% C.L.) |77j , 
i?(p 7 /K*7) = B(B -> fry)/B(B -> iT*7) < 0.28 (90% C.L.) TO] . 



(37) 
(38) 
(39) 



In the numerical analysis we will use only the last constraint since the ratio of branching 
ratios is theoretically cleaner as only the ratio of the form factors is involved, which is 
calculable more reliably. Concentrating on the neutral 5-decays, the LO expression for 
R(fn/K*i) is 



2B{B° -> p°(770) 7 ) 



where 



B(B -> K*^) 



V, 



td 



v t 



ts 



Ml- Ml* 



C s 7 {m b ) 



Ft 



(40) 



(41) 



with Fi ' K being the form factors involving the magnetic moment (short-distance) 
transition; M B = 5.2794 GeV, M K * = 0.8917 GeV, M p o = 0.7693 GeV and f = 0.58 [0. 

{m%) is the Wilson coefficient of the magnetic moment operator for the transition 
b — > (d, s) computed in the leading order approximation. Note that the annihilation contri- 
bution (which is estimated at about 25 % in B ± — > p ±r y |80|,|8l[|) is suppressed due to the 
unfavourable colour factor and the electric charge of the d-quark in B°, and ignored here. 
In the SM, these two Wilson coefficients coincide while, in the SUSY model we consider, 
they differ because of the effect of the insertion 5% i : 



C S 7 {M W 
C^{M ] 



w ) 



C 7 W {M W )+C 7 H (M w ) 
C^M W )+C^ /H (M V 



C 7 X (M, 



w ) 



(42) 
(43) 
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where , and C* are, respectively, the SM, the charged Higgs and the chargino 

contributions and their explicit expressions can be found for instance in Refs. P3L^]. The 
explicit expression for the mass insertion contribution is 



Cf J (M, 



w , 



K 



ad 



v t 



Id 



i=l 



6M3 

Xi 



COS 6*r 



X, 



i2 K^tiXt' x qxJ ( °u L t 2 



COS /?£ 

= °7 °«Lt 2 



(44) 
(45) 



where the loop functions ff n and f^ 11 are given in Appendix B. Using Eqs. ( |4*3] ) and ( fI5D 
it is possible to rewrite the ratio R(py/K* , y) in the following way: 



R{prt/K*i) = R 



SM 



le-p^MI 



U L t2 



Cf(m 6 ) 



(46) 



in which C 7 and Cj(rrib) are both real, 1723 is a QCD factor numerically equal to 0.66 and 
we have abbreviated R{fry / K*"f) SM by R SM for ease of writing. 



4 Numerical Analysis of SUSY Contributions 

In this section we study the correlations between the possible values of the parameters / and 
g as well as the numerical impact of b — > S7, 5 a + and b — > with the following experimental 
constraints 

2.41 x 10" 4 < B(B -> X s7 ) < 4.02 x 10~ 4 (95% C.L.) 

10 x 10~ 10 < 5a^+ < 74 x 10~ 10 (95% C.L.) (47) 

R(pj/K*j) < 0.28 (90% C.L.) 

where we have used the estimates of 5a /J +(SM) from Eq. (|30P . We perform the numerical 
analysis by means of high density scatter plots varying the SUSY input parameters over the 
following ranges: 



n = 


(100 4- 


1000) GeV , 


M 2 = 


(100 4- 


1000) GeV , 


tan (3s = 


3 4-35 




M H ± = 


(100 4- 


1000) GeV , 


M i2 = 


(100 4- 


600) GeV , 


h = 


-0.3 4 


-0.3 . 



(48) 



Notice that, according to the discussion of the previous section, we restrict the scatter plot 
to positive fi values only. Negative values are strongly disfavoured both by the 5a^+ > 
bound and by the b — > 37 branching ratio. In fact, it is possible to show that if \x is negative, 
the chargino contributions to b — > 57 tend to interfere constructively with the SM and the 
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charged Higgs ones. In order not to exceed the experimental upper limit, a quite heavy 
SUSY spectrum is thus required. In such a situation, high / and g values are quite unlikely. 

In Fig. ^|we plot the points in the (/, \g\) plane that survive the b — > 57, 5a^+ and b — > 
constraints. Scanning over the parameters given in Eq. (|4q) we find that the constraints in 
Eq. (f|7]) restrict / and \g\ to lie essentially in the range / < 0.4, \g\ < 2.0 . We also find 
that the sign of C x /CY tt is positive over all the SUSY parameter space that we scanned. 

The impact of 5a^+ on our analysis is not very strong, once we limit the scanning to the 
/i > region only. Moreover, as follows from Eq. (^), the size of the chargino contribution 
to 5a^+ is controlled by the mass of the muon sneutrino. In our framework, mp is a free 
parameter and we have imposed the lower bound 5a^+ > 10 x 10~ 10 in the loosest possible way 
by choosing = 100 GeV (a value that is reasonably safe against direct search constraints). 
On the other hand, we can not reject points which give a too large contribution to 5a M + 
because a large enough sneutrino mass can always suppress the SUSY diagram and reduce 
5a^+ to a value smaller than 74 x 10~ 10 . Notice that, if > 300 GeV, all the points that we 
consider do satisfy the upper bound: in order to obtain larger contributions it is necessary 
to impose a very light sneutrino mass. Only models in which the squark and the slepton 
mass spectra depend on the same inputs will be able to fully exploit the correlation between 
the anomalous magnetic moment of the muon and observables related to B physics. 

The impact of the b — > d*y constraint is taken into account by imposing the following 
upper bound on the mass insertion: 



°U L h 



< 1 ^iS- 1 



C?(m 6 ) 



r]23C 7 



(49) 



Again, we find that with the current experimental bound R{p r )/K* , y) < 0.28, most of the 
otherwise allowed (/, \g\) region survives. This situation will change with improved limits 
(or measurements of R(pry/K*j)). 

In Fig. |3] we perform the same analysis presented in Fig. ||] but we allow only for points 
that give a positive sign for the Wilson coefficient Cj computed in the LO approximation. 
The issue whether it is possible or not to change the sign of Cj depends on the model 
and has been long debated in the literature. In particular, this sign strongly characterizes 
the behaviour of the forward-backward asymmetry and of the dilepton invariant mass in 
b — > s£ + £~ transitions, as well as the sign of the isospin violating ratio A (see below) and of 
the CP-violating asymmetry in the radiative decays B — > fry |79j . We use Cf for calculating 



the CP-asymmetry in B — > fry. The quantity A in the NLO approximation requires the 
Wilson coefficient C^ LO . However, as shown in Ref. |[79|| , A is stable against NLO vertex 
corrections. Recently, also the so-called hard spectator corrections have been calculated to 
0(a s ) in B — ► p7 [0,|83[] with the result that A is stable also against these corrections. We 
refer to Refs. |5^ , [79| , |84| -|86[] for a comprehensive review of the positive C-j phenomenology. In 



Fig. ^ open circles represent points that satisfy the b — > 37 and 5a^+ constraints. The black 
dots show what happens when the experimental bound R(fry / K*^) < 0.28 is imposed. In 
implementing this constraint we use Eq. ([49]). If for a given point 5 1[m turns out to be smaller 
than 1, we plot |gf|^_^im, otherwise we set 8 = 1. It is important to note that the dependence 
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of Cj H and C^ 11 on the mass insertion is, respectively, linear and quadratic. From Fig. |3] 
one sees that all the points that are compatible with a positive C-j provide, indeed, a too 
large contribution to b — > dj, and hence are effectively removed by the cut on R(pj / K*j). 
This result is quite reasonable because, in order to change the sign of CV, a large positive 
chargino contribution is needed: since C* and C 7 depend on the same input parameters we 
expect their magnitude to be closely correlated. In Fig. |], we show explicitly the correlation 
between C 7 {ra^) and C 7 (m b ) = r]2sC 7 (M w ) in both the negative and positive Cf allowed 
regions. In the second plot, in particular, \C 1 M1 '{m^)\ turns out to be greater than one for all 
the points: this implies that the mass insertion constraint is always non trivial. The strong 
bound shown in Fig. |3] is obtained by taking into account that g depends quadratically on 
5u L t 2 - The conclusion is that if C 7 > is experimentally established, our analysis implies 
strong constraints on the quantity \g\ from b —>■ dj decays, permitting only small deviations 
from the MFV- value: \g\ = 0. 



5 Unitarity Triangle Analysis 

In section ^ we have shown that the impact of this class of SUSY models on observables re- 
lated to the unitarity triangle (UT) can be parameterized by two real parameters and by one 
phase (see Eqs. (H) and (|10D ). In this section we analyze the implications of this parametriza- 



tion on the standard analysis of the UT. As usual we use the Wolfenstein parametrization |36 
of the CKM matrix in terms of A, A, p and r\: 

( 1-f A AX 3 {p-t7])- 

A 1-f AX 2 | . (50) 



V 



\A\ 3 (1- p-irj) -AX 2 1 

In the following analysis we extend this parametrization beyond the leading order in A; 
as a consequence it is necessary to study the unitarity triangle in the plane (p, ff) where 
p = p(\- X 2 /2) and fj = r] (1 - A 2 /2) 0. 

Let us collect the relevant formulae for ex, A^B (ds) and a^x s as functions of /, g and 
5 c L t 2 - 

e K = - G \^^ K M K ^ lm {^VccC^ + 2KX; Vtc Cr c + K\tCf«{l + f + g)}e%, 

(51) 

Gj 
6n' 2 

a^K s =sin2(/3 + d ) (54) 

where X q = V qd V* s , = |arg(l + f + g), (3 denotes the phase of V t * d and from Eq. (|50| ) it 
follows 

i - p 2 ) + ri 2 
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AM B d = -^VBM B J 2 Bd B Bd M 2 w \V t ^; d \ 2 CY tt \l + f + g\, (52) 
AM Bs = -^ B M B J 2 B B B M 2 w \V t ^:fCY tt {l + f) , (53) 



Parameter 



Value 



Vtt 

Vcc 
The 

B K 
Vb 



f. 



B B d 



0.57 88 



1.38 ±0.53 
0.47 ±0.04 



0.94 ±0.15 f91],[92 



0.55 88 



230 ± 40 MeV 



i.i6 ±0.05 g 



Table 1: Parameters used in the UT-fits 

The quantities %, r] cc , r] cc , and are NLO QCD corrections. Their values together with 
those of the other parameters are collected in Table |I[ 

Our first step is to investigate the regions of the parameter space spanned by /, and 
gi that are allowed by the present experimental data. The procedure consists in writing 
the x 2 °f the selected observables and in accepting only values of / and g which satisfy the 
condition min p r) (x 2 ) < 2. In the computation of x 2 we use the following input 

e K = (2.271 ±0.017) 10~ 3 [ PDG [11]], (56) 
AM Bd = 0.484 ± 0.010 ps" 1 [HFWG |T|], (57) 

= 0.090 ±0.025 [PDG 0], (58) 



V, 



ub 



v cb \ 

a ^ Ks = 0.79 ±0.12. (59) 

The value quoted above for AM Bd is the present world average of experiments on B - B 
mixings, including the recent BABAR || and BELLE || measurements. We introduce the 
experimental data on AM Bg in calculating the chi-squared (x 2 ) using the so-called amplitude 
method f95] . The prescription consists in adding to the x 2 the term x 2 am b = (1 — A) 2 l a \ 
where A is the amplitude of the (B s — B s ) oscillation, given by (1 ± ^4 cos AM Bg t), and 04 
is the corresponding error. Both A and 04 are functions of AM Bs - Notice that using this 
method, the statistical interpretation of the value of the x 2 m its minimum is preserved. In 
Ref. the authors include the AM Bg data using an alternative procedure. They consider a 



log-likelihood function referenced to AM Bs = 00 and add the term A log L°° = (1/2 — A) /a a 
to the x 2 ■ m this way the significance of the data is increased. Notice that, in order to 
interpret the output of this method in terms of confidence levels it is necessary to perform 
a monte-carlo based analysis [p7 |. Since we are interested in the statistical meaning of the 



minimum of the x 2 1 we prefer to use the standard amplitude analysis. 

We present the output of this analysis in Fig. |5|a. For each contour we fix the value of 
/ and we require the xL to be less than 2. Since, as can be seen in Fig. ^ / is always 
smaller than 0.4, we restrict the analysis to / = 0, 0.2 and 0.4. Moreover, for each of these 
values, we require \g\ not to exceed the upper limit which, according to the analysis of the 
SUSY contributions presented in Fig. 0, we set respectively to 1, 2 and 1.5. The equation 
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sin2(/3 + 9 d ) = a^Ks nas two solutions (mod n) in which (/? + 6 d ) lies in the ranges [0, 7r/4] 
and [vr/4, tt/2] respectively. This implies that, given a value of /, we expect two distinct 
allowed regions in the (gn, gi) plane, that we call respectively small and large angle solutions, 
characterized by f3 + 9 d < vr/4 (i.e. 6 d ~ O(0. 1) and [3 + 6 d > vr/4 (i.e. 6 d ~ 0(1)). These two 
regions have some overlap since in the limit a^,K s = 1 (which is allowed at 90% CL.) the two 
solutions coincide. We find that, once the upper bound on \g\ is imposed, only a small part 
of the / = 0.2 and a tiny corner of the / = 0.4 large angle solution survive. Improvements 
in the experimental determination of the B — > X s 7 and B — > fry branching ratios as well 
as more stringent lower bounds on the SUSY spectrum will have a strong impact on the 
allowed \g\ values. On the other hand, more precise measurements of (V^/V^l, a^ s , AM Bs 
and progress in the determination of the relevant hadronic parameters, will contribute to 
reduce sizably the size of the allowed (gx, gi) regions. In view of these considerations, we 
expect the large angle (8 d ~ 0(1)) solution is less likely to survive in future and we will 
concentrate in the following on the small angle scenario, i.e., 9 d ~ O(0.1). 

It is interesting to note that the amplitude method, that is conservatively used to set the 
constraint AM Bs > 14.9 ps _1 at 90% C.L., also yields a 2.5 a signal for oscillations around 
AMb s = 17.7 ps _1 . This would-be measurement is equivalent to a determination of the 
ratio AMb s / AMs d which in turn depends on the precisely computed hadronic parameter £ s : 
its impact on the unitarity triangle is thus expected to be quite significant fl5fl . In Fig. |b 



we assume this signal to be a measurement with AMb s = 17.7 ± 1.4 ps _1 in order to explore 
its implications on the previous analysis. The effect of the assumed value of AMb s is to 
reduce the gn, gi > allowed regions; moreover, its impact is stronger the higher is the value 
of /. Note, in particular, that, with the assumed value of AMb 3 , the / = 0.4 contour almost 
disappeared. This happens because the experimentally favoured high values of a^ Ks tend to 
sharpen the mismatch between the AMg d and AM^/AM^ constraints (which is due to a 
non-zero value of /). 

Before concluding this section we would like to show the impact of the Extended-MFV 
model considered in this paper on the profile of the unitarity triangle in the (p, fj) plane, 
and the corresponding profiles in the SM and MFV models. In Fig. || the solid contour 
corresponds to the SM 95% C.L., the dashed one to a typical MFV case (/ = 0.4, g — 0) and 
the dotted-dashed one to an allowed point in the Extended-MFV model (/ = 0, g% = —0.2 
and gi = 0.2). The representative point that we consider survives all the experimental 
constraints examined in the previous section. Using the values of p and fj that correspond 
to the central value of the fit, we obtain the following results for the various observables: 
\V ub /V cb \ = 0.079, e K = 2.27 x 10" 3 , AM Bd = 0.484 ps" 1 , AM Bs = 20.7 ps^ 1 and a^ Ks = 
0.81. If \g\ is sufficiently large, 9 d can be regarded as an essentially free parameter and the 
fit will choose the value that gives the best agreement with the experimental measurement. 
In Fig. |7| we plot the CP asymmetry function of arg<5 Sit ~ 2 (expressed in degrees). 

The light and dark shaded bands correspond, respectively, to the SM and the experimental 
1 a allowed regions. The solid line is drawn for / = and \g\ = 0.28. The experimental 
band favours arg^^ in the range [0°, 100°]. Employing the explicit dependence 

9 d = 1 arg(l + / + \g\e 2iw ^Lh) ( mo d vr) , (60) 
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Contour 


9R 9i 




AM Bs 


a^Ks 


a 7 


1 


0.2 0.2 


0.094 


20 ps- 1 


0.78 


119° 40° 


2 


0.0 -0.2 


0.110 


20 ps" 1 


0.71 


101° 51° 


3 


-0.4 0.1 


0.081 


17 ps" 1 


0.73 


64° 98° 



Table 2: Central values of the CKM ratio \V u b/V c b\, the B s — B s mass difference AMb s , the 
CP asymmetry a^K s and the inner angles a and 7 of the unitarity triangle for the contours 
in the extended-MFV model plotted in Fig. [|. The {gn,gi) values for the contours are also 
indicated. 

the above phase interval is translated into —3° < 6d < 8°, for the assumed values oi\g\ and 
f, which is a typical range for 9d for the small angle solution with the current values of a^x s - 

In order to illustrate the possible different impact of this parametrization on the unitarity 
triangle analysis, we focus on the / = case and choose three extremal points inside the 
allowed region in the plane (gn, gj). We concentrate on the small angle scenario. In Fig. | we 
plot the 95% C.L. contours in the (p, fj) plane that correspond to the points we explicitely 
show in Fig. ^|a. We summarize in Table || the central values of \V u b/V c b\, AMb s , the CP 
asymmetry a^K s , the inner angles a and 7 of the unitarity triangle computed for the different 
contours. Contour 1 is drawn for positive gj and 9d is consequently positive. This implies 
that a^K s is expected to be larger than in the SM. The CP asymmetry and IV^/V^I are 
very close to their world averages while the angle 7 is smaller than •y SM . In contour 2 
the phase 9d is negative and the CP asymmetry is thus predicted to be lower than the 
experimental central value and still 7 < r y SM . Contour 3 is drawn for g = —0.4 + O.li and 
is particularly interesting since it corresponds to a solution in which a^K s is larger than in 
the SM and the Wolfenstein parameter p is negative, i.e. p < 0, implying a value of the 
inner angle 7 in the domain it/ 2 < 7 < it. This is in contrast with the SM-based analyses 
which currently yield 7 < 7r/2 at 2 standard deviations [T2"|-|H| and with the other solutions 



shown in Fig. [8|. We note that analyses p8| , |49[| of the measured two-body non-leptonic 
decays B — > tttt and B — > Kir have a tendency to yield a value of 7 which lies in the range 
7 > 7r/2 (restricting to the solutions with fj > 0). While present data, and more importantly 
the non-perturbative uncertainties in the underlying theoretical framework do not allow to 
draw quantitative conclusions at present, this may change in future. In case experimental 
and theoretical progress in exclusive decays force a value of 7 in the domain 7r/2 < 7 < n, 
the extended-MFV model discussed here would be greatly constrained and assume the role 
of a viable candidate to the SM. 

In Fig. [5] and Table |3] we show the consequences, on the analysis described above, of 
reducing the error on the CKM ratio IK^/V^I by a factor of 2. The value IV^/V^I = 
0.090±0.013, used by us to illustrate the improved constraint on the CKM unitarity triangle, 
is obtained from the inclusive measurement of \V c b\ from CLEO JO] and LEP ||98|| , yielding 



\Vcb\ = (40.57 ± 1.21) x 10 , and the average of the currently measured values of \V u b\ by 
the CLEO |U and LEP [(§] groups, yielding \V ub \ = (3.64 ± 0.46) x 10~ 3 . Note that the 
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Contour 


9R 9i 




AM Bs 


a^Ks 


a 7 


1 


0.2 0.2 


0.092 


20 ps- 1 


0.77 


120° 40° 


2 


0.0 -0.2 


0.102 


20 ps" 1 


0.65 


101° 53° 


3 


-0.4 0.1 


0.084 


17 ps" 1 


0.75 


64° 97° 



Table 3: The same as in Table^. The error on \V u b/V c b\ is reduced by a factor of 2. 



reduced error on IK^/V^I does not affect sizably the existence of the p < solution in the 
extended-MFV model shown here. 



6 Implications of the Extended-MFV model for b —> dj 
transitions 

In this section we study the implications of the extended-MFV model on the observables 
related to the exclusive decays B — > p 7 , namely the ratio R(p r ) / 'K*"f) defined in Eq. (|40|), 



the isospin breaking ratio 



where 



and the CP asymmetry 



A± o = ^42)_ l! (62) 



2B(B° -> p° 7 ) 



A = B(B~ -> p - 7 ) ~ B(B+ -> p+ 7 ) 
CP B(B- -> p - 7) + -> p+ 7 ) ' 1 ' 

We perform the numerical analysis for following set of SUSY input parameters (that satisfy 
all the constraints previously discussed): p = 120 GeV, M 2 = 350 GeV, tan/? = 4, M^ 2 = 
280 GeV, 6i = —0.29 and M H ± = 290 GeV. This allows us to exploit in detail the dependence 
of the various observables on the phase of the mass insertion and to give an illustrative 
example of the modifications in the profile of these quantities. 



In Fig. [T^ we plot the ratio R(p^ / K*^) as a function of aig5^ 2ilL in the Extended-MFV 
model, and compare the resulting estimates with the SM estimates, shown by the two dashed 
lines representing the la SM predictions. The solid curves are the SUSY results obtained 
for p and rj set to their central values and for |5„ L t 2 | = <5 lim , 0.55 and 0.25. Here, 5 lim is the 
absolute value of the mass insertion that saturates the experimental upper bound R < 0.28; 
it is required to be smaller than 1 and it depends on the phase of the mass insertion. For 
the point that we consider, it varies between 0.6 (for arg<% 2 £ L = 0,n) and 1. The shaded 
region shown for the |5g f 2 | = <5 lim represents the 1 a uncertainty in the CKM-parameters 
(p, fj) resulting from the fit of the unitarity triangle. In the maximal insertion case, the 
experimental upper bound R(pj / K*j) < 0.28 is saturated for arg(^ 2Si G [0, 7r/2]U[37r/2, 2tt}. 
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Note that if we require the absolute value of the insertion to be maximal, the ratio R{p r y /K* 1 -)) 
is always larger than in the SM. We point out that, in the extended-MFV model, this ratio 
does not show a strong dependence on the p and i] values as long as these CKM parameters 
remain reasonably close to their allowed region. On the other hand, the impact of reducing 



\5$ % \ is quite significant. 



Taking into account the discussion at the end of Sec. [|and that Cj{m b ) is negative for all 
the points that allow for a sizable mass insertion contribution, the arg5j 2fi region in which 
the experimental bound on R(p n f / K*j) is saturated turns out to be strongly dependent on 



the sign of C^ 1 . 



Moreover, a large C 



MI 



is usually associated with a large stop mixing angle 

In 



whose sign determines, therefore, the overall sign of the mass insertion contribution 
our case, C 7 < and the region Re5^ 2ilL > is consequently favoured. The qualitative 
behaviour of this plot can be understood rewriting Eq. fl46|) as 



R{p~f/K*~f) = R SM \1 + Ae 4arg5 | 2 = R SM [1 + A 2 + 2v4cos(arg<5)] 



(64) 



where A = IS^^r/^^^C^ 11 /C 7 (m b ) is positive for the point we consider. 

The explicit expressions for the isospin breaking ratio and the CP asymmetry in the SM 
79] 



arc 



A C p 



2e, 



Fx + -e A (F? + Ft) 



2F 2 (AJ - e A A\' 
Cf M K)(l + A LO ) 



(65) 
(66) 



where e A 



-0.3, A u j 



0.046, A\ 



(i)t 



-0.016 and 



F 2 



Re 



Im 



^ub^ud 

v^yFd 



ub v ud 



v tb v* d 

^ub ^ud 



v tb v* d 



cos Ot 



sin a 



(67) 
(68) 



Note that e A is proportional to l/Cf M (mh). Eqs. fl6q) and (|66f) can be easily extended to 
the supersymmetric case by means of the following prescriptions: 



Kd Kd ex p{* ar gC' 

■ \C^m b )\ . 



[i.e.: a — > a — argC^Wb)] 



C 7 SM (m b 



(69) 
(70) 



In Figs. and ^ we present the results of the analysis for the isospin breaking ratio and 
of the CP asymmetry in B ± —>■ p ± 7, respectively, for the three representative cases: \$u L i 2 \ = 
<5 lim , 0.55 and 0.25 in the Extended-MFV model and compare them with their corresponding 
SM-estimates. Since, in all likelihood, the measurement of the ratio R(pj / will precede 
the measurement of either A or the CP-asymmetry in B — > pj decays, the experimental value 
of this ratio and the CP-asymmetry a^ Ks can be used to put bounds on \5\ and arg^ j 2 . 
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The measurements of A and the CP-asymmetry in — > p ±r y will then provide consistency 
check of this model. Concerning the cases |^f 2 | = 0.55 and 0.25, we must underline that 
large deviations occur in the phase range in which an unobservably small R(pj / K*^) is 
predicted. On the other hand, it is interesting to note that, for the case of maximal insertion 
and for a phase compatible with the measurements of a^,K s (see Fig. [?[), sizable deviations 
from the SM can occur for the CP asymmetry but not for the isospin breaking ratio. 

7 Summary and Conclusions 

The measured CP asymmetry a^ Ks in B decays is in good agreement with the SM prediction. 
It is therefore very likely that this CP asymmetry is dominated by SM effects. Yet, the last 
word on this consistency will be spoken only after more precise measurements of a^K s and 
other CP- violating quantities are at hand. It is possible that a consistent description of CP 
asymmetries in 5-decays may eventually require an additional CP-violating phase. With 
that in mind, we have investigated an extension of the so-called Minimal Flavour Violating 
version of the MSSM, and its possible implications on some aspects of B physics. The non- 
CKM structure in this Extended-MFV model reflects the two non-diagonal mass insertions 
from the squark sector which influence the FCNC transitions b — > d and b — > s (see Eq. (pTT]) ) . 
In the analysis presented here, we have assumed that the main effect of the mass insertions 
in the 5-meson sector is contained in the b — > d transition. This is plausible based on the 
CKM pattern of the b — *> d and b — > s transitions in the SM. The former, being suppressed 
in the SM, is more vulnerable to beyond-the-SM effects. This assumption is also supported 
by the observation that the SM contribution in b — > 57 decays almost saturates the present 
experimental measurements. We remark that the assumption of neglecting the mass insertion 
(Sc L t 2 ) can be tested in the CP- asymmetry in the b — > s sector, such as Aqp(B — > X s j), 
Acp(B — > K*j), the B® - B® mass difference AMg s , and more importantly through the 
induced CP-asymmetry in the decay B®(B®) — > J/ip(f), which could become measurable in 
LHC experiments |37j due to the complex phase of (S £ j 2 ). The parameters of the model 
discussed here are thus the common mass of the heavy squarks and gluino (Mq), the mass of 
the lightest stop (M^ 2 ), the stop mixing angle (%), the ratio of the two Higgs vevs (tan 
the two parameters of the chargino mass matrix (fj, and M 2 ), the charged Higgs mass (Mh±) 
and the complex insertion (S^ i 2 ). 

We have shown that, as far as the analysis of the unitarity triangle is concerned, it is 
possible to encode all these SUSY effects in the present model in terms of two real parameters 
(/ and gn) and an additional phase emerging from the imaginary part of g (gi)- We find 
that despite the inflation of supersymmetric parameters from one (/ in the MFV models) 
to three (in the Extended-MFV models), the underlying parameter space can be effectively 
constrained and the model remains predictive. We have worked out the allowed region in the 
plane (/, \g\ = \gn + igi\) by means of a high statistic scatter plot scanning the underlying 
supersymmetric parameter space, where the allowed parametric values are given in Eq. fl48|). 
The experimental constraints on the parameters emerging from the branching ratios of the 
decays B — > X s ^ and B — > /ry (implemented via the ratio R(pj/K*j)), as well as from the 
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recent improved determination of the magnetic moment of the muon {g — 2) M were taken 
into account, whereby the last constraint is used only in determining the sign of the //-term. 

We have done a comparative study of the SM, the MFV-models and the Extended-MFV 
model by performing a x 2 -analysis of the unitarity triangle in which we have included the 
current world average of the CP asymmetry a^x s (Eq. @) and the current lower bound 
AMg s > 14.9 ps _1 using the amplitude method. Requiring the minimum of the x 2 to be 
less than two, we were able to define allowed-regions in the gi) plane (correlated with 
the value of /), which are significantly more restrictive than the otherwise allowed ranges 
for \g\. We studied the dependence of the CP-asymmetry a^ Ks on the phase of the mass 
insertion (S% i 2 ) and find that, depending on this phase, it is possible to get both SM/MFV- 
like solutions for a^K S i as well as higher values for the CP-asymmetry. We constrain this 
phase to lie in the range 0° < arg5 Si ^ 2 < 100°, which typically yields the Extended-MFV 
phase 9d to lie in the range —3° < 6 a < 8°. The assumed measurement of the mass difference 
AMb s , when inserted in the \ 2 analysis, further restricts the allowed regions in \g\. However, 
as AMg s has not yet been measured, this part of the analysis is mostly illustrative. Finally, 
we have shown the profile of the resulting CKM-Unitarity triangle for some representative 
values in the extended-MFV model. They admit solutions for which a^ Ks > and 
7 > 7 SM , favoured by present data. 

To test our model, we have focused on three observables sensitive to the mass insertion 
(Su L i 2 ) related to the radiative decays B — ► pj. We have worked out the consequences of 
the present model for the quantities R(pj)/R(K*j) = 2B(B° -> p°j)/B{B° -> K*°j), the 
isospin violating ratio A ±0 = B(B ± — > p ±r y)/2B(B° — > p 0r y) — l, and direct CP-asymmetry in 
the decay rates for B — > p7 and its charge conjugate. We conclude that the partial branching 
ratios in B — > pj, and hence also the ratio R(p , y/K*'-f) can be substantially enhanced in this 
model compared to their SM-based values. The CP-asymmetry can likewise be enhanced 
compared to the SM value, and more importantly, it has an opposite sign for most part of 
the parameter space. On the other hand, it is quite difficult to obtain a significant isospin 
breaking ratio without suppressing the branching ratios themselves. 

Finally, we remark that our analysis has led us to an interesting observation: the require- 
ment of a positive magnetic moment Wilson coefficient (i.e., Cy > 0), entering in b — > 57 and 
b — > s£ + £~ decays, is found to be incompatible with a sizable contribution to the parameter 
g, which encodes, in the present model, the non-CKM flavour changing contribution. Thus, 
it is possible to distinguish two different scenarios depending on the sign of C|. In the C| < 
case, as in the SM, only small deviations in the b —>■ s phenomenology are expected, but 
sizable contributions to gn and gi are admissible, thereby leading to striking effects in the 
b —>■ d sector. On the other hand, a positive C7 will have strong effects in the b — > s sector 
but, since \g\ will be highly constrained, the model being studied becomes a limiting case of 
the MFV models. In particular, in this scenario, no appreciable change in the CP-asymmetry 
0"ipK a compared to the SM/MFV cases is anticipated. Since the experimental value of a^K s 
[Eq. (g)] is in agreement with the SM/MFV-models, configurations in which a^,K s receives 
small corrections have a slight preference over the others. These aspects will be decisively 
tested in S-factory experiments. 
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A Stop and chargino mass matrices 

The 2x2 stop mass matrix is given by 



2 / M~ M~ \ 

* \M? M~ ) 1 ] 

v tLR tRR. ' 



where 



Mf LL = M| + (1 - ^ sin 2 6 W ) cos 2(3 S Mf + M 2 , (72) 

Mf RR = M~ + ^ sin 2 6 W cos 2(3s M§ + M 2 , (73) 

= M t |A-//*cot(/3 s )| . (74) 

The eigenvalues are given by 



2M h 2 = ( M L + M D ± \/( M L- M L) 2 + 4 ( M L) 2 > ( 75 ) 

with M 2 2 < M~ . We parametrize the mixing matrix 1Z 1 so that 

tl\ rrtil h \ ( cos0 { sm9 { \ ( t L 



t 2 J 11 { t R ) [ -an0 { cosOi ) { t R ) ' (76) 
The chargino mass matrix 

M x + _ ( M 2 M w V2smp s \ 

~ [ M w V2cos(3 s p ) {77) 

can be diagonalized by the bi-unitary transformation 

U* a M$V k * p = M n 5 jk , (78) 

where U and V are unitary matrices such that M~+ are positive and M~+ < M~+. 

X,j Ai A. 2 
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B Loop functions 



The loop functions for box diagrams, entering in ek, AMr, and AMb b , are, 



G(a,b) 
Yi(a, b, c, d) 



ab ( a 2 -8a + 4 , b 2 - 8b + 4 , , 

in a + — — 7T7T in o — 



Y 2 (a, b, c, d) 



4 \(a-b)(a-l) 2 

(6 — a)(c — a)(d — a) 
| c 2 _ 

(a — c)(b — c)(d — c) 
a 



(b-a)(b-l) 2 

b 2 

{a - b){d - b){d - b) 
d 2 



(a-l)(6-l). 



x 4cd 



(a — d)(b — d)(c — d) 
6 



In 6 
In d 



(b — a)(c — a)(d — a) 
c 



In a + 



lnc + 



(o-6)(c-6)(d-6) 
d 



ln6 



ir"(a,&, c, d) 



(a — c) (6 — c) (d — c) (a — d) (6 — d) (c — d) 

Yi(a, a, c, d) + Yi(b, b, c, d) — 2Y\(a, b, c, d) 



In d 



(a - 6) 2 

The loop functions for penguin diagrams, entering in b — > (s, d)7 and in the anomalous 
magnetic moment of the muon, are 

x 3 — 6a; 2 + 3a; + 2 + 6x In x 



Fx(x) 
F 3 (x) 

h{x) 



\2{x - l) 4 
x 2 - Ax + 3 + 2 In x 

2(x- l) 3 ' 
-7 + 12a; + 3a; 2 - 8x 3 + 6x(-2 + 3x) log a; 

6(x- l) 4 
5 - 12x + 7a; 2 - 2z(-2 + 3x) log a; 

2(x- l) 3 ' 

A,2(a;) - /i,2(y) 
(a; - y) 
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Figure 1: Feynman diagrams contributing to C* and C± respectively. The bubble represents 



the mass insertion 8, 
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Figure 2: Allowed points in the (/, |g|) plane. The ranges of the supersymmetric parameters 



and the constraints from b — ► 57, £0^+ and 6 — > c/7 t/iat we impose are given in Eq. ( \4Bi ) and 
Eq. respectively. 
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Figure 3: Allowed points in the (/, \g\) plane that are compatible with a positive value of 
the Wilson coefficient Cj. The empty circles satisfy the b — ► 37 and 5a^+ constraints. The 
impact of imposing in addition the upper bound on the b —>■ transition is represented by 
the black dots. 




Figure 4: Correlation between the SUSY contributions to Cjim^) and C 7 (nib). The plots 
correspond to the negative (left figure) and positive (right figure) solutions for C% (nib) allowed 
by the experimental bounds on B(B — ► X s j). Note that Cj(m b ) = Cj(m b ) + Cj H (nib) ^u L t 2 - 
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Figure 5: (a) Region of the (gn,gi) plane for which mm p ^(x 2 ) < 2 /rom t/ie CKM-UT 
fits. The AMb s constraint is taken into account using the amplitude method. The contours 
correspond to f = 0, 0.2 and 0.4 and the constraints on \g\, coming from Fig. ^, are \g\ < 1, 
2 and 1.5 respectively. The shaded areas correspond to a solution with 9 + (3 > 45° and 
are allowed for f = 0.2 and f = 0.4. The three dots are representative points that we use 
to illustrate the impact of the parametrization we propose on the unitarity triangle, (b) 
The same as in (a), but interpreting the 2.5 o enhancement in the amplitude as AMb s = 
17.7 ± 1.4 ps-\ 
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Figure 6: Allowed 95 % C.L. contours in the (p,fj) plane. The solid contour corresponds to 
the SM case, the dashed contour to the Minimal Flavour Violation case with (f = 0.4, g — 0) 
and the dashed-dotted contour to the Extended- MFV model discussed in the text (/ = 0, Qr = 
-0.2, # = 0.2). 





Figure 7: The CP asymmetry a^K s as a function o/arg5 s £ 2 expressed in degrees. The solid 
curve corresponds to the Extended- MFV model (/ = 0, \g\ = 0.28). The light and dark shaded 
bands correspond, respectively, to the allowed 1 o region in the SM (0.58 < a^K s < 0.82,) 
and the current 1 o experimental band (0.67 < a$K s — 0.91/ The plot on the right shows 
the correlation between arg^^ and the angle 6d = |arg(l + / + |fi , |e 2 * arg s £*a) j (modn). 
The experimentally allowed region favours 0° < arg<5 SLt ~ 2 < 100° that translates into —3° < 

e d <8°. 
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Figure 8: Allowed 95% CL. contours in the (p,fj) plane. The solid contour is the SM case. 
The two semicircles represent the 2 o region allowed by \V u b/V c b\ = 0.090 ± 0.025. The 
contours numbered 1 to 3 correspond respectively to the points indicated in Fig. [5]a, and their 
{QRiQi) values are given in Table |M Note that the UT-contour 3 yields 7 > tt/2. 




Figure 9: The same as in Fig. [|. The error on \V u b/V c b\ is reduced by a factor of 2. The 
two semicircles represent the 2 a region allowed by \V u b/V c b\ = 0.090 ± 0.0125. The (gn,gi) 
values corresponding to the three extended-MFV contours are given in Table [|. 
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Figure 10: The ratio R(py /K*^) = B(B — ► py)/B(B — > K*j) as a function of arg <5j 2fii 
(in degrees) in the Extended-MFV model, satisfying the present experimental upper bound 
i?(p7 / 'K*j) < 0.28 (at 90% C.L.). The solid lines are obtained for p andfj set to their central 
values and for f | = 5 Um , 0.55 and 0.25. The shaded region in the top curve represents 
the 1 a uncertainty due to the fit of the unitarity triangle. The dashed lines indicate the 1 o 
SM prediction. 
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Figure 11: The isospin breaking ratio A as a function of arg<5f 2Si (in degrees). See the 
caption in Fig. WR for further explanations. 
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Figure 12: The CP asymmetry in — > p ± 7 as a function of axg(5t 2ilL ) (in degrees), 
the caption in Fig. [77] for further explanations. 
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